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Experimentally observable Quantum Accelerator Modes are used as a test case for the study of 
some general aspects of quantum decay from classical stable islands immersed in a chaotic sea. The 
modes are shown to correspond to metastable states, analogous to the Wannier-Stark resonances. 
Different regimes of tunneling, marked by different quantitative dependence of the lifetimes on 1/h, 

r .i. are identified, depending on the resolution of KAM substructures that is achieved on the scale of fi. 

The theory of Resonance Assisted Tunneling introduced by Brodier, Schlagheck, and Ullmo is 

^— ^ revisited, and found to well describe decay whenever applicable. 
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0\ ■ I. INTRODUCTION. 

' Classical Hamiltonian systems generically display mixed phase spaces, where regions of regular motion and regions 
ÇNJ , of chaotic motion coexist 0. Although no classical transport is allowed between different, disjoint regions, quantum 
C~-» ■ transport is made possible by so-called " dynamical tunneling" , which allows wave-packets to leak through classical 
invariant curves. In the case of a wave-packet initially localized inside a classical stable island, a transfer of probability 
into the chaotic region arises, that may continue a long time, and thus take the form of irreversible decay from the 
island, whenever h is so small that the Heisenberg time related to fast chaotic diffusion is much larger than the period(s) 
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, of regular motion in the island. Islands of regular motion have been found to play crucial roles in many contexts, 
such as optical cavities Q, driven cold àtoms 0, where they have inspired detailed theoretical analysis P|, and 
billiards , where they motivated mathematical investigations [í§ and atom-optical experimental realizations [J^ Wide 
^p^ 1 attention has been attracted by Chaos Assisted Tunneling, which denotes transmission between symmetry-related 
"tí | islands, across a chaotic region Q. Slow tunneling out of islands dominates the dynamical localization properties in 
^ ■ some extended systems @, • 

In this paper we explore decay from regular islands into the surrounding chaotic regions, which is a central theoretical 
CJ^ 1 issue in all the above hinted subjeets. It is an appealing idea that decay rates may be estimated, using only information 
drawn from the structure of classical islands. On account of the ubiquitous character of mixed phase spaces, this 
theoretical programme is attracting significant attention 0, 0, lïïl UJ. In this paper we address the problem in a 
special case, which has a direct experimental relevance. The Quantum Accelerator Modes (QAM) were experimentally 
discovered when cold caesium àtoms, falling under the action of gravity, were periodically pulsed in time by a standing 
wave of light [Ï3L IÏ5| . The theory of this phenomenon shows that the dynamics of an atom is described, in an 
appropriate gauge, by a formal quantization 28] of either of the classical maps : 

Jt+i = Jt + ~ksm(0 t+1 )±2nn , t +i = 6 t ± J t mod(27r) . (1) 

The classical phase portraits exhibit periòdic (in J) chains of regular islands. The rest of phase space is chaotic, 
and QAMs are produccd whenever an atòmic wavepacket is at least partially trapped inside the islands. Therefore, 
QAMs eventually decay in time, due to quantum tunneling [ï^. The problem of QAMs has a relation to the famous 
Wannier-Stark problem |l6| about motion of a particle under the combined action of a constant and a periòdic in 
space force field 0], an d the class of QAMs we consider in this paper are due to metastable states 0, which are 
analogous to the Wannier-Stark resonances, and are associated with sub-unitary eigenvalues of the Floquet evolution 
operator. We derive and numerically support quasi-classical, order-of-magnitude estimates for their decay rates, 
based on classical phase-space structures. We theoretically and numerically demonstrate that decay is determined by 
different types of tunneling, depending on how significant the KAM structures inside the island are on the scale of 
H, and in particular we show that these different mechanisms result in different quantitative dependence on the bàsic 
quasi-classical parameter given by A/h, where A is the area of an island. In the regime where higher-order structures 
inside the island are quantally resolved, we use the theory of Resonance Assisted Tunneling [íj to describe the decay 
of QAMs and find a good agreement in the presence of a single dominant resonance. In particular we observe, in 
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especially clean form , a remarkable stepwise dependence on the quasi-classical parameter, as first predicted in [91. 
The phase space islands which are studied in this paper are directly related to experiments on làser cooled àtoms [Ï3l 
lïfij . Unfortunately the step structure which is here predicted occurs in parameter ranges where the decay rate is 
extremely small. Finding parameter ranges where decay rates are appreciable, and still exhibit step dependence, is a 
challcnging task for experimental application. 



II. CLASSICAL AND QUANTUM DYNAMICS. 



QAMs one-to-one correspond to stable periòdic orbits of the map on the 2-torus which is obtained from (IJ on reading 
J modulo (2w) |líj|l4(. This class of orbits in particular includes fixed points (i.e., period-1 orbits) of map 0J. It 
is these period-1 orbits that give rise to the most clearly observable QAMs, and in this paper we restrict to them 
|29). Map differs from the Standard Map only because of the drift 2tïVL in the lst equation. Despite its formal 
simplicity, this variant introduces nontrivial problems, concerning Hamiltonian formulation and quantization, which 
are discussed in this section. 



A. Wannier-Stark pendulum. 

Let k = ke and 2-kVI — ae, with e a small parameter. For e = 0, map (JIJ has circles of fixed points at the "resonant" 
vàlues of the action J = 2i:s, (s € Z). Straightforward calculation shows that for e > a stable fixed point, surrounded 
by a stable island, survives near each resonant action, whenever k is larger than a and smaller than a stability border 
|l4j . Other islands , related to periòdic orbits of higher periods, may or may not significantly coexist with such 
pcriod-1 islands, depending on parameter vàlues. In any case, numcrical simulation shows that motion outside all 
such islands is essentially chaotic at any e > 0. 

Using canonical perturbation theory, one finds [l6| that in the vicinity of resonant actions, and at lst order in e, the 
dynamics Q are canonically conjugate to the dynamics which are locally ruled by the "resonant Hamiltonian": 



H res (J,ti) = -J 2 +eV(i?) ; V{&) = + fccos(i?) . 



(2) 



Multivaluedness of V($) is removed on taking derivatives, so Hamilton's equations are well-defined on the cylinder, 
even though the Hamiltonian J2J is not. They uniquely define a "locally Hamiltonian" flow on the cylinder, which will 
be termed the Wannier-Stark (WS) pendulum, because, if d were a linear coordinate and not an angle , then J3J) would 
be the Wannier-Stark (classical) Hamiltonian for 1-dimensional motion of a partició in a sinusoidal potential combined 
with a static elèctric field |Ï7j. Trajectòries of the Wannier-Stark pendulum are obtained, by winding around the 
circle the trajectòries which are defined on the line by the Wannier-Stark Hamiltonian. 

If a = 0, @ is the Hamiltonian of a Standard pendulum, and motion is completely integrable in either of the two 
regions in which phase space is divided by the pendulum separatrix. If a ^ this is not true anymore. In particular, 
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FIG. 1: Stable and unstable manifolds of a Wannier-Stark pendulum (|5J with ae = 0.5 and ke = 0.7. On the cylinder, arrows 
lie on a single continuous trajectory. 
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if \a\ < \k\ then the flow has one stable and one unstable fixed point. Motion is completely integrable inside the stable 
island dclimited by the separatrix, but not outside, because trajectòries are unbounded there and motion cannot be 
confined to a torus (Fig^}. The separatrix (also called a "bounce" 0] or "instanton trajectory" [23|) is the trajectory 
which approaches the unstable point, both in the infinitely far past, and in the infinitely far future ( loop in Fig^ linc 
PQ in FigEJ ■ The resonant Hamiltonian provides but a local description of the motion near a resonance. It misses 
the periodicity in action space which is an important global feature of the problem; nevertheless, it does provide a 
description of the inner structure of stable island(s). 

B. Quantization. 

Quantization of map Q is a nontrivial task, because a shift in momentum by 2nil, as in the lst eqn.Q, may be 
inconsistent with quantization of momentum in múltiples of h. This problem disappears if the angle 8 in £3> is replaced 
by x g] — oo,+oo[, because then the map describes motion of a particlc in a line, and straightforward quantization 
yields the unitary operator 

1j e iatX / h e -ik<i cos(X) / h e -iP 2 /2h 

where X and P are the canonical position and momentum operators. However, the quantum dynamics thus defined 
on the line do not define any dynamics on the circle, except in cases when U commutes with spatial translations by 
2tt. This case only occurs when ae is an integer múltiple of fi. Then quasi- momentum is conserved and Standard Bloch 
theory yields a family of well-defined rotor evolutions, parametrized by vàlues of the quasi-momentum. If ae = mh/n 
with m and n intcgcrs, then it is casy to sec that hi n commutes with spatial translations by 2tt and so the n— th power 
of the classical map "on the circle" can be safely quantized. 

Similar subtleties stand in the way of quantizing the Wannier-Stark pendulum. The WS Hamiltonian "on the line" 
never commutes with translations by 2-7T, as long as ae ^ 0. However, if ae = mh/n, then the unitary evolution 
generated by the WS Hamiltonian over the integer time n does commute with such translations |18| . and so it yields 
a family of unitary rotor evolutions. Each of these yields a quantization of the WS pendulum flow at such integer 
times. We shall restrict to such " commensurate" cases. In the language of the theory of Bloch oscillations these 
are the cases when the "Bloch period" Tb = h/(ae) and the kicking period are commensurate. 

III. DECAY RATES. 

Let U generically denote the unitary operators that are obtained by quantization of map Q or powers thereof, as 
discussed in Sect III Bl We contend that, despite classical stable islands, the spectrum of U is purely continuous. Quasi- 
modes related to classical tori in the regular islands correspond to metastable states, associated with eigenvalues of Ü, 
which lie strictly inside the unit circle, and thus have positive decay rates T . They are analogous to the Wannier-Stark 
resonances. Arguments supporting this contention are presented in Appendix^ along with methods of numerically 
computing decay rates T. In this Section we obtain order-of-magnitude estimates of decay rates T. 
Motion inside the islands is not integrable, but just quasi-integrable, and displays typical KAM structures, such as 
chains of higher-order resonant islands. We separately consider the cases when such structures are small (resp., large) 
on the scale of h. In the latter case, our bàsic theoretical tool is the notion of Resonance Assisted Tunneling, which 
was introduced in This theory is presented from scratch in sects. IIII Bl and IIII Cl with special attention to the 
role of classical and quantum perturbation theory. 

A. Wannier-Stark tunneling 

First we consider the case when h is small compared to the size of an island and yet large compared to the size of 
the stochastic layer and of resonant chains inside the island. This in particular means that 2nd order corrections 
on the resonant Hamiltonian @ are classically small, and so one expeets the bare resonant Hamiltonian to capturo 
the essential features. This Hamiltonian is formally the Wannier-Stark Hamiltonian and stable islands are associated 
with elliptic motion near the bottom of the potential wells (Fig|2l . A WKB estimate for the smallest decay rate from 
a well is 



r _ ^0 e -2S(A,B)/h 
27T 



(3) 
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where luo is the angular frequency of the small oscillations and iS(A, B) is the imaginary action along the classically 
forbiddcn path from point A to point B, at constant energy equal to the value of the potential at the bottom of the 
well: 

s(a,b) = f B de VHv(e) - v{e A )) . 

Je A 

Rcficction 6 — + tt — 6 turns 9a hito 9q, 9b into 9p, and reverses the sign of the argument of the square root; and so 
S(A, B) = S(P, Q), the real action along the path from P to Q at constant energy equal to the value of the potential 
in Q. This path is the separatrix, so 2S(A, B) is equal to the area enclosed by the separatrix, which is in turn nearly 
equal to the area A of the actual island in the regime we are considering. Therefore, 

T ~ f£ e -*/* (4) 

This result is compared with a numerical simulation in FigEI Even better agreement with numerical data is obtained 
by using in 10 the trajectory with energy hui /2 above the bottom of the well, which is an approximation to the 
ground state energy in the harmònic approximation. This is shown in Figs^and^l 

The success of the elementary WKB approximation is due to the fact that in cases like Figs[3]and2]thc potential 
barriers on the right of -ds are significantly lower than V(i!}a), so that tunncling trajectòries have to cross just one 
potential barrier. When this condition is not satisfied, one is faced with the full complexity of the WS problem. 



B. Phenomenological Quantum Hamiltonian. 

In Fig0{b) we show the dependence of T vs fi, for the case of the island of FigQJa). Like in the casc of Fig|3 2nd 
order resonances are not quite pronounced here, and so, at relatively large vàlues of fi (leftmost part of the Figure) 
the dominant contribution to decay is given by WS tunneling, and good agreement is observed with the theory of 
sect IIII Al As fi decreases, a clear crossover is observed to slower decrease of T, indicating that a different mechanism 
of decay is coming into play which overrules WS tunneling. In our interpretation, this is the mechanism of "resonance 
assisted tunneling" to be discussed in what follows, which was introduced in ref. 9J. It is a quantum manifestation 
of classical KAM structures, and so we now consider the case when fi is small compared to the size of the island, 
yet phase-space structures produced by higher-order corrections to the resonant Hamiltonian are not small on the 
scale of fi. The approach to be presently described is based on a quantum Hamiltonian, which is not formally derived 
from the exact dynamics, but is instead tailored after the actual structure of the classical phase space. This heuristic 
approach is applicable when the WS linewidth (0} is negligible with respect to the coupling between different WS 
tori (and between WS tori and the continuum) which is due to higher order corrections, and in fact it totally misses 
WS tunneling given by It is assumed that the classical partition island/chaotic sea is quantally mirrored by a 
splitting of the Hilbert space of the system in a "regular" and a "chaotic" subspace, with respective projectors P r 
and P c ; and that the Hamiltonian H may be written as 

È = H r + H c + V + í> f , (5) 

where H r = P r HP r is a "regular" Hamiltonian, H c — P C HP C is a "chaotic" Hamiltonian, and V = P r HP c couples 
regular and chaotic states. In the case of maps, the Hamiltonian formalism has to be recovered by means of Floquet 
theory. Eigenvalues of Floquet Hamiltonians fall in Floquet zones, and for a system driven with a period T in time, 
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FIG. 3: (Color online).(a) Phase portrait of system Q with negative sign, for k = 0.8 and 2iríl — 0.7. (b) Circles: numerically 
computed decay rate from the center of the regular island, as a function of 1/h. Lines show WKB estimates, obtained by using 
classically forbidden paths at constant energy @. These are equal either to a minimum of the WS potential (dashed line) 
shown in Fig|21 or to the ground state energy in the potential well, estimated in the harmònic approximation (sòlid line). 



the width of a zone is 2nh/T. In our case we assume T = 1 and so we identify the lst Floquet zone with the interval 
[0, 2-kK). For the Hamiltonian H c we assume in our case (where a ^ is always understood) a continuous spectrum 
|30| As this Hamiltonian is assumed to be " chaotic" , we further assume that its quasi-energies are non-degenerate 
and completely fill each Floquet zone, because typical random matrices drawn from circular ensembles have simple 
spectra, with eigenphases uniformly distributed in [0, 2ir}. The discrete eigenvalues E n of the regular Hamiltonian H r 
are therefore immersed in the continuous spectrum of the chaotic Hamiltonian H Cl so the coupling perturbation V 
drives them off the real axis, and the imaginary parts —ihT n /2 they acquire are estimated by Fermi's Golden Rule : 

T„ ~ ^\(E,c\VÏ\E n )\% =En . (6) 

where \E, c) is the eigenvector of H c associated with an eigenvalue E in the continuous spectrum (5— function normal- 
ization in energy is assumed for such eigenvectors). Choosing E n in the lst Floquet zone of H c , © may be rewritten 
as 

r„ ~^||p vt|^)||V(^n). 

n 

where Pq is projection onto the lst zone of H c . The function /j,(E) is the Local Density of States, normalized to 1, of 
the vector PoV^\E n ) with respect to the chaotic Hamiltonian H c . It yields the probability that a transition prompted 
by V> from state \E n ) will lead to the continuum state \E,c) in the lst zone. One may introduce the "ergodic" 
assumption, that all in-zone transitions have the same probability. Then 

r„ ~ L·hv^Ejw*. (7) 
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FIG. 4: (Color online).Same as FigEUfor k = 0.7 and 2vrfi = 0.5. The island chain highlighted in (a) is generated by a 11 : 1 
resonance at Zm « 0.251, with u = 4.376 • 10 -8 , and M = 2.785 (see App. [HJ. The dotted line in (b) represents the theory 
l|16fl for decay assisted by this resonance, using the unperturbed spectrum l|10|l . 

Whether or not the ergodic assumption is accepted, (0) may be assumed to hold up to a factor of order 1. 



C. Resonance-assisted decay. 

At lst order in e, the Hamiltonian H r should correspond to the classical resonant Hamiltonian, and its eigenstates to 
quantized tori thereof. Therefore, the coupling V only reflects classical corrections of higher order than lst, because 
the resonant Hamiltonian has no coupling between the inside and the outside of an island. Higher order corrections 
are present in H r as well, and adding their secular (averaged) parts to the resonant Hamiltonian a new Hamiltonian 
H roB is obtained, which is still completely integrable, and shares the action variable X of the integrable WS pendulum 
flow (inside the WS separatrix). Semiclassical quantization of H IQS yields energy levels e„ = H IBS (T n ) with quantized 
actions X n = (n + l/2)h. 

In the following, by " perturbation" we mean what is left of higher-order corrections, after removing averages. Thus 
the "unperturbed" quasi-energy eigenvalues of H r are given by e n + 2irNH with JV any integer ( 2ttH is the width 
of a Floquet zone). In the classical case, the destabilizing effects of the perturbation are mainly due to nonlinear 
resonances. A classical r : s (r, s integers) resonance oceurs at X = X rs if roj(I rs ) = 2tts, where co(X) = dH Ioa (X)/dX 
is the angular frequeney of the integrable motion in the island. The quantum fmgerprints of classical nonlinear 
resonances are quasi-degeneracies in the quasi-energy spectrum. A degeneracy appears in the unperturbed quasi- 
energy spectrum of H ri anytime two or more energy levels e„, e m , ... are spaced by múltiples of 2irh, that is the width 
of a Floquet zone. In the vicinity of a classical resonance X rs , this is approximately true, whenever n — m is an integer 
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múltiple of r; for, in fact, |e„ — e m | ~ h\n — m\u>ÇZ rs ) — 2Trhs\n — m\/r. Thus, in the semiclassical regime, a classical 
resonance induces quantum quasi-degeneracies, which involve whole sequences of quantized tori. As such tori are 
strongly coupled by the perturbation, decay is enhanced. This quantum effect is Resonance Assisted Decay Q. 
The perturbative approach assumes that the sought for metastable states basically consist of superpositions of such 
strongly coupled, quasi-degenerate states. This leads to replacing P r in (0 by projection P R onto a quasi-degenerate 
subspace, similar to methods used in [27|: and so, in order to use Fermi's rule (0, P R HP R must be diagonalized, and 
P R V must be specified. To this end we first write the matrix of P R HP R in a basis of quasi-degenerate states. Let T no 
be a quantized action close to X rs . The energy levels e ri0 +Nr, where the integer N takes both negative and positive 
vàlues, have an approximately constant spacing, close to 2nhs, and thus form a ladder. There is one such ladder for 
each choice of the integer no in the set of the r closest integers to n ra = T rs /h — 1/2 (not necessarily an integer), and 
we fix one of them. The quasi-energy levels e„ 0+ jvr — 2nhNs are quasi-degenerate . Denoting \N) the corresponding 
eigenstates, the projector P R onto the ladder subspace is a finite sum J2n l-^X-^li where N ranges between a N* < 
and a N* > 0. The numbcr of terms in the sum is equal to the number of nearly resonant states inside the island 
and so is given by L + 1 where L — N* — N* . Using that the actions Jjv of levels in a r : s ladder are approximately 
spaced by múltiples of rh, L is estimated by 



L « Int 



\2nh j r 



(8) 



If only nearest-neighbor transitions are considered, the matrix of P R HP R over the basis \N) is tridiagonal, of size L+l. 
The off-diagonal elements v(N) — (N\H\N + 1) may be semiclassically assumed to slowly change with N, and will be 
hence denoted simply by v. The diagonal elements are the nearly degenerate quasi-energies W(N) — e no 

+Nr- Nsh, 

and Taylor expansion of ií ros (1) to 2nd order near T r8 yields 

— s h 2 

W(N) » P rcB (X rs ) + 2i:h-(n - n rs ) + — (n + rN - n rs ) 2 (9) 

r ZM 

where M — l/u)'(I rs ). It follows that, apart from a constant (independent on N) shift, the diagonal matrix elements 
of P R HP R are 

W(N)^^(rN + Sn) 2 . (10) 

where Sn = uq — n rs . Replacing this in P R HP R , one easily recognizes that the classical limit h — > 0, L — > oo of P R HP R 
is the classical pendulum Hamiltonian: 

H rs (I*,e*) = r ^ r +2vcos(9*) , (11) 

in appropriate canonical variables X*,8*(Nh ~ I* = (2" —T rs )/r as h — > 0). This is the well-known pendulum 
approximation near a classical resonance [lj , directly derived from quantum dynamics|22j , and is related by a simple 
canonical transformation to the slightly different Hamiltonian I|B1J) which is used in 

The coupling to continuum P R V remains to be specified. Of all quantized tori in the chain, the closest to the chaotic 
sea corresponds to state \N*), and we assume that this one state (in the given chain) is coupled to the continuum. 
This assumption implies : 

VÏP R = |x> (AH , (12) 

where \x) is some vector in the chaotic subspace, about which our one assumption is that it lies in the lst Floquet 
zone. Its norm has the meaning of a hopping amplitude from the "gateway state" \N*} to the normalized state 
The latter state may be thought of as a "last beyond the last" nearly resonant state, corresponding to 
an unperturbed torus which was sunk into the stochastic sea by the perturbation and so cannot any more support 
a regular quasi-mode of H r . Thus one may denote \\x)\\ = v(N*), and extrapolate to this last transition, too, the 
semiclassical assumption v(N*) ~ v. After all such additional constructions, the "ergodic assumption" (end of section 
IIH B|l is just that this last state has equal projections on all eigenstates of H c . Fermi's rule Q now yields, for the 
decay from an eigenstate \E m ) (0 < k < L) of P R HP R l 

r m ~ l ^\(N*\E m )\ 2 . (13) 
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The labeling m = 0, . . . , L of the eigenstates of P R HP R is arbitrary for the time being. The tridiagonal Hamiltonian 
is defined on a chain of states and the scalar product in the last formula is the value of the eigenfunction (N\E m ) 
at the rightmost site N = N* in the chain. Assuming the eigenfunction to attain its maximum modulus at a site 
N m , its value at site N* should be of order exp(— Ç m \N* — N m \). The quantity Ç^ 1 is the fall-off distance of this 
eigenfunction. For a tridiagonal Hamiltonian on a chain, arguments by Herbert, Jones, and Thoulcss |2l| . estimate 
this distance , as 

U - In (^) , (14) 



where D m is the geomètric average of the differences \Ej — E m \ (m fixed, j variable, j ^= m), that is 

1 L 

\n(D m ) = - H\E,-E m \), (15) 



and v is the geomètric average of the hopping coemcients v(N), (JV* < N < N* — 1). Under the assumption 
v(N) ~ v =const., v — v, and so finally 



2 / \ 2(N*-N m ) 

v í v x 



It should be emphasized that the E m used in 115|l are the eigenvalues of P R HP R and not its diagonal elements W m 
(fTÜ|> . and that eqn. ((TU) is not perturbative. 



D. Single-Resonance Assisted Decay. 

At any v ^ 0, an island hosts a dense set of resonances, but only a few of them are quantally resolved. On the 
other hand, quasi-resonant ladders of states may be formed, only if the resonant transitions are not so broad as to 
involve off-ladder states. That means that the classical chain of islands should not be too wide, because its width is 
determined by the same parameter v which yields the hopping amplitude between nearly resonant states. 
In this section we consider the case when a single, not too wide resonance dominates all the others, so that each 
metastable state may be assumed to sit on one of the ladders which are associated with that resonance. Decay rates 
r m are not afhxed to unperturbed tori, but to eigenstates \E m ) of a ladder Hamiltonian. For v — these correspond 
to quantized tori in the quasi-resonant ladder, via (|1U|) . and their labeling by m — 0, . . . , L may be chosen accordingly. 
For not quite small v, however, the correspondence between metastable states and unperturbed tori may be broken, 
due to avoided crossings; and so labeling by the original quantum number may not anymore reflect how deep in the 
island an eigenstate is located. 

We restrict the following discussion to decay " from the center of the island" , meaning that we consider a metastable 
state, which is mostly supported in the innermost part of the island, and is labeled by m = 0. For this state, l|ltj|) 
reads: 

„,2(L + 1) 2 

r ° ~ -Wïïf = ¥ e ~ 2 '° L 5 & = MA.) - (17) 

For quite small v, explicit calculation is possible, using for Ej their unperturbed (v = 0) vàlues given by (|10(l . This 
leads to 

\E Ó - £ | w h 2 \r 2 j 2 - 2n rs rj\/2M, (18) 
and then l(Tü|l with STty and (jT5|> leads to 

e ° - ln (2^) + 7 è ln d^' 2 - 2n ^n\) ■ (19) 

v / j=1 
to be used in l|17|) . In the average, (|17|l takes the form (cp.JSJl): 



r ~ v e -^ A '* rn , . (20) 
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This equation shows that the average dependence of Tq on í/h is exponential, because the dependence of Ço on h is 
semiclassically weak. In fact, v is a classical quantity, and in the limit S->0 the average over levels Ej which enters 
eqn. lfTÏÏ)l turns into a purely classical quantity, given by the phase-space average of ln(Tí rs — E ) (cp.JH}) over the 
island. Approximating the sum in (|19|l by an integral, and denoting A rs = 2iïT rs = 2T:hn rs the area enclosed by the 
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FIG. 5: (Color online).(a) Phase portrait of map Q with negative sign, for k = 2.5 and 2iríl — 1. The 4-islands chain is 
due to a 4 : 1 resonance at T41 ~ 0.43, with v ~ 7.275. 10 -4 and M « 3.866. (b) Minimal decay rate from the island shown 
in (a), calculated by different numerical (symbols) and theoretical (lines) methods. Circles: truncated basis method. Pluses: 
long-time decay of an initial coherent state in the center of the island. Small Squares : Complex Scaling Method. Sòlid line : 
formula 1171 . with 1191 . Dotted-dashed line: continuum approximation 12ÜL12H . 

r : s resonant unperturbed torus, and x — Ars/·A, 

Ço ~ - ln ( !f V ^\ - 2 + 2x\n(2x) + (1 - 2a) Infli - 2x|) . (21) 



V A 2 

In formula (|17|l (with (|19|l ). L is anyway a discrete variable, which discontinuously jumps by 1 anytime increase of 
í/h grants accommodation of a new quasi-resonant torus in the island (ep.©)- This produces a stepwise dependence 
of Tq on í/h, superimposed on the average exponential dependence. This structure is smoothed in (|2Ü|I with l|21|) . 
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E. Numerical results for resonance assisted decay. 

Eans. H17|) and (|20|l (supplemented by (|19fl or JUJ) are the main results of the theory of resonance assisted decay and 
are tested against numerical results in this section. They require specification of M, v,T rs as input parameters. In the 
case of well-pronounced chains of resonant islands, the vàlues of M, v, and I ra may be found by "measuring" areas 
in the classical phase-space portrait, taking advantage of the fact that a resonant chain is bounded in between the 
separatrices of a pendulum (|llf> . This method was introduced in 0, and is reviewed in Appendix IBI for the reader's 
convenience. 

For small v, unperturbed eigenvalues IjlOII may be used in (|15fl for the purpose of calculating D m , leading to (|17|) (with 
(JEU), but one has to beware of vàlues of h, that enforco degeneracy of the unperturbed spectrum 1101) . Inadvertent 
use of the spectrum (jlCÍ|> in (|15|) in such cases leads to artificial divergence of T. In the actual spectrum to be used 
in Ijl5(l . degeneracies are replaced by avoided crossings, which may lead to local enhancement of resonance-assisted 
tunneling, as discussed in Appendix [Ç] and shown in Fig. 1111 

A case with a single dominant resonance 4 : 1 is presented in FigO . Decay rates calculated from l|17f) with l|10|) and 
(|19fl are shown in (b) for the innermost state in the island which is shown in Fig|SJa). Also shown are results of 
directly calculating Ts, by methods described in Appendix 1X1 Formula l|17fl . with H19|) . is seen to correctly reproduce 
the actual Ts, in order of magnitude at least. The stepwise dependence predicted in and explained in the end of 
previous section is here remar kably evident. The continuum (quasiclassical) approximation (|20|) . H21|) is also shown in 
the same Figure. In Fig|H|it is seen to better and better agree with l|17(l . with Ijl9|) . at smaller and smaller vàlues of H. 
As remarked in section IÏII Dl the presented theory, being essentially perturbative, is expected to fail in the case of 
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FIG. 6: (Color online).Comparison between EJ, with l|19^ . and the continuum formula )fi . (121 fl 

large classical resonances. An example is presented in Figd Finally in FiglHlwe present a case with two resonant 
chains of comparable size. Results are not well described by the theory based on either resonance, and thus appear to 
contradict a somewhat natural expectation, that each resonance should contribute its own set of metastable states. 
The single-ladder picture may not be adequate in such cases, which therefore remain outside the present scope of the 
theory. 

IV. DISCUSSION AND CONCLUSIONS 

The decay rates of some metastable states related to phase space islands were calculated and the required theoretical 
framework was developed. The main results of the paper are (gj, (|Po)l and (|T7|) . where the decay rates of wave packets 
in phase space islands were calculated for various conditions. If the effective Planck's constant is sufficiently large, 
so that island chains cannot be resolved on its scale, Standard WKB theory was found to work well. As the effective 
Planck's constant is decreased, island chains are resolved, and dominate the decay by the mechanism of resonance 
assisted tunneling. Its signature here is the step structure of Fig[5fb). The average slope of lnT (see l|20[l ) as a 
function of í/h is — ^A/nr, with £o given by (|21|l and is independent of h. In the WKB regime a different slope 
is found (see (@J). For some vàlues of the parameters, resonance-assisted tunneling can be further enhanced by a 
degeneracy between a semiclassical state deep inside the island and one that is close to the boundary. An interesting 
question is about possible effects of this kind, due to "vague tori" l2q. i.e.. to classical structures which quantally act 
as tori, in spite of lying outside an island. Indeed, observations in [ÏJ| have suggested a possible role for cantori in 
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FIG. 7: (Color online) .Illustrating a case with a very broad resonant chain. (a): Phase portrait of MapQ with negative sign, for 
k = 7r and 2-kQ, = 0.5. The numerical vàlues for the 3 : 1 resonance are /31 = 0.13, v = 1.9 • 10 -3 , M = 1.52 (b): Minimal decay 
rate from the island, vs 2iï/h. Circles: numerical data (truncated basis). Dotted, dashed, and dash-dotted lines show decay 
rates for the three lowest-lying states in the island, computed from the theory of decay assisted by the large 3 : 1 resonance . 
The sòlid line shows the smallest of these. 

enhancing QAMs at times. 

The theory strongly relies on the dominance of one resonant island chain If the phase space area of the chain is not 
small as is the case in Fig|7]or if the tunneling is assisted by two (or more) resonant island chains of approximately 
equal strength, as is the case in FigO °ur theory requires modification. 

The theory is relevant for systems of experimental interest but the steps of FigJS] were found for a regime where 
the decay rate is too small to be experimentally accessible. Overcoming this problem is a great theoretical and 
experimental challenge. 
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APPENDIX A: METASTABLE STATES. 

Our methods of computing decay rates T were: (1) basis truncation, (2) complex scaling, (3) simulation of wave-packet 
dynamics. In the cases investigated in this paper, the most economical one, and thus the one of our prevalent use, 
was (1). The other two methods were used to cross-check results of (1) in a number of cases. The observed agreemcnt 
between such completely independent computational schemes demonstrates the existence of resonances (in the sense 
of metastable states). 
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FIG. 8: (Color online) .Illustrating a case with two classical resonances of comparable size. Here k = 1.329 and 2iríl = 0.5336. 
Circles in (b) represent the minimal decay rates from the truncated basis method. Lines in (b) represent the theory of Resonance 
Assisted decay, using either the 6 : 1 resonance (sòlid line) or the 7 : 1 resonance (dashed line). The former resonance has 
v = 1.8 • 10~ 4 , M = 4.504, lei = 0.46, and the latter has v = 3.1 ■ 10~ 4 , M = 2.626, 171 = 0.93. 



1. Basis Truncation. 

Let U denote the unitary evolution operator that is obtained by quantizing map ijljl. as described in sect IIIBl and 
let \n), (n 6 Z) be the eigenvectors of the angular momentum operator J, such that (9\n) = (27r) -1 / 2 exp(in&). 
"Basis truncation" consists in replacing U by U v = P^UPv, where P v = Y^\ n \<v \ n )( n \- This introduces an artificial 
dissipation, which turns the quantum dynamics from unitary to sub- unitary. The eigenvalues Zj, (j = 1, . . . , 1v + 1) 
of U v lie inside the unit circle, with positive decay rates Tj = ln(l/|zj|)/2. As v is increased, most of them move 
towards the unit circle, but some appear to stabilize at fixed locations inside the circle, because they approximatc 
actual, subunitary eigenvalues of the exact (y — 00) non-dissipative dynamics. The seeming contradiction to unitarity 
of the limit dynamics U is solved by the observation that the eigenfunctions of U v , which are associated with such 
eigenvalues, tend to increase in the negative momentum direction (Fig llO(l . as expected of Gamov states, and if this 
behavior is extrapolated to the limit, then they cannot belong in the Hilbert space wherein U acts unitarily. In 
order to make room for such non-unimodular eigenvalues, U must be extended to a larger functional space. Complex 
sealing, to be described in the next subsection, provides a consistent method of doing that. 

Simulations of wave packet dynamics, performed in the total absence of any dissipation whatsoever, confirm this 
interpretation of the stable subunitary eigenvalues (PigJSJ) . The initial wave packet is a coherent state supported near 
the center of the island . We use a FFT algorithm, so the computed evolution is fully unitary, and reliably reproduces 
the exact evolution over a long time, thanks to the large dimension of the FFT . We compute the decay in time of the 
probability in a momentum window which contains the classical island. After an initial ràpid decay, due to escape 
of the part of the distribution which initially lies outside the island, the decay turns into a clean exponential with 
rate 7 w 5.66 ÍFig llOfl . Among the eigenfunctions of the truncated basis evolution, which correspond to stabilized 
eigenvalues, we select the one which has largest overlap with the chosen initial state. We thus find that (i) the decay 
rate T of this eigenfunction is rs 7, (ii) the Husimi function of that part of the wavepacket, which has survived in the 
chosen window until the end of the dynamical calculation, nearly reproduces the Husimi function of the eigenfunction 
(Fig0. 

Numerical computation of long-time exponential decay may be less easy than in the above particular example. If an 
initial state is overlapped by several metastable states with slightly different T, resolving them may take quite a long 
computational time, and hence a huge basis, because of accelerated motion outside the island. 
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FIG. 9: (Color online).(a),(b),(c): Evolution of an initial coherent state located at the center of the island of Fig|^] for k = 2.5, 
2-kQ, = 1, h = 0.25, numerically simulated by a FFT method with basis size 2 17 . (a),(b),(c) show contour plots of the Husimi 
functions at times t=100, t=1000, t=16000 respectively. (d): Husimi contour plot of the eigenfunction of the truncated basis 
dynamics, which has the largest overlap with the initial coherent state. The size of the truncated basis is 4096. 

2. Complex Scaling. 

Scattering resonances may be sometimes computed by diagonalizing a non-hermitean Hamiltonian, which is con- 
structed by "complex coordinate" methods such as analytic dilation, and the like [Ul^j]. Despite absence of scat- 
tering theory, a method of this sort was devised for the subunitary eigenvalues considered in this paper, as follows. 
Any function ip(ff) over [0, 2n] is at once a function ip{ z ) of the complex variable z running on the unit circle. Let 
B denote the class of those functions ip(0), which can be analytically continued to the whole complex plane, except 
possibly the origin. For given 1 > p > the scaling operator h p is defined to act on the functions of this class as 
in (h p ip)(9) = i/j(pe ). The crucial property of the operator U, which makes the present construction possible, is 
that of transforming functions in class B in functions in the same class, and so the operator U p = h p Uh' p l is a well 
defined operator in B. This operator trivially extends by continuity to an operator U p |25j | which is defined on the 
whole of L 2 ([0, 2tt]). This extension formally amounts to defining U also on a class of functions, which are not square 
integrable. The new "functions" thus acquired in the domain of the evolution operator may be very singular objects; 
e.g., in the momentum representation, they are allowed to exponcntially diverge at infinity. In the special case h=l/n 
with n integer, U p has the form : 

Ü p = (0 » e i«e e -^+™ s (é) e fc-si„(è) e - l J 2 /2 ) k± = ik\p±p~ 1 \ . 

which restitutes U for p = l. The eigenvalues of U p as a L 2 operator are at once eigenvalues of the "extended" U, and 
each of the latter eigenvalues is an eigenvalue of U p for p sufficiently distant from 1. The "complex scaling" method 
of computing T, which is mcntioncd in the main text, consists in diagonalization of U p . At small H this method is 
computationally problemàtic, due to exponcntially largc clements in the matrix of U p . 

APPENDIX B: INFERRING PARAMETERS OF A CLASSICAL RESONANCE FROM PHASE 

PORTRAITS. 

Parameters T rs , M, and v of a classical r : s resonance respectively specify the value of the unperturbed action wherc 
the resonance is located, the inverse nonlinearity l/ui'(T rs ), and the strength of the resonant harmònic perturbation. 
These parameters are indispensable for the formalism described in sect llll Cl and may be retrieved from the phase- 
space portrait, using formulae taken from [Tï| . Here we reproduce a sketchy derivation for the reader's convenience. 
Motion in a r : s resonant chain is approximately described by a pendulum Hamiltonian [lj , which may be written in 
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FIG. 10: (Color online).(a): Decay in time of the probability in a momentum window containing the island, for the same 
parameter vàlues as in Fig|ïï] and for the same choice of the initial state. The size of the basis used in the simulation is 2 18 . 
The straight line has slope « —5.66. (b) squared modulus of the truncated basis eigenfunction presented in Fig[ÜJd), in the 
momentum representation. 




7.2 7.4 7.6 7.8 8 8.2 8.4 8.6 8.8 9 -2 -1 



FIG. 11: (Color online).(a) Decay rate from the lst excited state in the island, vs 1/h, for the same parameter vàlues as in FiglSI 
Dots: numerical (truncated basis method). Sòlid line: formula I16H with 1151 1. using the unperturbed ladder spectrum I1UI . 
Squares: same formula, using the perturbed spectrum of the tridiagonal matrix with I1U1 as diagonal elements. Dashed line: 
same formula, using the actual spectrum, semiclassically reconstructed from the layout of tori in the island. (b) the eigenstate 
of the (truncated) evolution operator, corresponding to the decay rate shown in (a) for vàlues of 1/h: 7, 8, 8.2 and 8.6 (from 
above down), in the momentum J representation. Only a J— interval of roughly the size of the island is shown. (c) Decay rate 
r (circles) and real quasi-energy w (squares) vs 1/h, for the lst (full symbols) and the 5th (empty symbols) excited states. 



the form (which is canonically equivalent to 1)110: 

H rs {l,y) = _(X-Z„) 2 + 2v cos(r^) . (Bl) 

In this paper, 1 is the action variable of the WS-pendulum Hamiltonian J2J) . The Separatrices of Hamiltonian I|B1|) 
are the curves X = I±{tp) — T rs ± 2y / Mv(l — cos(r<^)), and so the phase areas S± they enclose satisfy: 



S+ + S- = 4ttX 



rs ; 
2tt 



S+ - S- =2 dcp (I+(<p) -T rs ) = 16V2Mv . (B2) 



The monodromy matrix of the stable period-r orbit that is responsible for the resonant chain is easily obtained 
by linearizing the flow I|B1|) near the stable equilibrium point(s). Its trace is found to be Ai = 2cos(ri^) where 
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v = r^/2v/M is the angular frequency of the small pendulum oscillations. This leads to 



yj2v/M = r- 2 arccos(Al/2) . (B3) 

The phase-space areas S± and the Monodromy Matrix can be numerically determined, and once their vàlues are 
known cas. (|B2l) and l|B3|) can be solved for l rs , M, and v. 



APPENDIX C: AVOIDED CROSSINGS. 



An exact degeneracy arises in the unperturbed ladder spectrum 11U|) whcnever two quantized actions in the ladder 
are symmetrically located with respect to the the resonant action X rs . This requires Sn = or Sn = ±r/2 in 11U|I and 
so, if (rii + í/2)H is the smallest quantized action in the ladder, such symmetric pairs exist if, and only if, 

\ = 2nt tZ +1 ^ m ~ 1 ' (C1) 

for some integer l > 1. In the above inequality, the length L of the ladder depends on h as in (JHJ). In the case of 
FigOO where = (the ground state), I|C1|I is never satisfied for í/h < 20. The data in Fig|S]were computed using 
the spectrum (|lü|) . and no significant difference could be found between them and data computed by using the actual 
spectrum, semiclassically reconstructed from the layout of tori in the island. In Fig llll where rij = 1 (the " lst excited 
state"), <|C1|1 is satisfied for í/h = 8.2, 12.8, 17.44, corresponding to / = 1,2,3. At such vàlues of í/h the spectrum 
p()[l(with the appropriate Sn) is degenerate and using it in formulas . obviously causes T to diverge, as shown 
by the narrow peaks in Fig^] Such artifacts disappear on inserting the proper (perturbed) spectrum, obtained by 
diagonalizing the ladder Hamiltonian with the diagonal elements given by (|10f) . because for d / the spectrum is 
never degenerate. Nevertheless avoided crossings take place at the vàlues l|Cl|l of h, giving rise to local peaks in the 
dependence of T on h. 

At the same vàlues of í/h avoided crossings are observed even between subunitary, stabilized eigenvalues of the 
truncated evolution operators (see App[X|. This is shown in Fig ll·lf ch The eigenvalues are written z = e ~ r / 2 ~™ 
and it is seen that, as í/h approaches a value ta 8.2, a pair of complex eigenvalues undergo a close avoided crossing. 
The corresponding states exhibit Standard behavior at avoided crossings. The distribution in momentum J of one of 
them is shown in Fig lll[) . This state nominally corresponds to the n = 1 unperturbed state, and in fact in (b)(top) 
it looks similar to the lst excited state of a harmònic oscillator. The other state nominally corresponds to the n = 5 
unperturbed state. At the avoided crossing (3d inset from top in (b)) the former state significantly expands over the 
island, because it is basically a superposition of two unperturbed states, which are located symmetrically with respect 
to 241 ■ This gives rise to a local enhancement of the decay rate. 
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